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^ ; The formula UmAr_^oo 'Ep<N,p,p+2 both prime log(p) log(p+2) = C2 is tested on the 

C^ . computer up to A^ = 2^^ ~ 1.1 x 10^^ and very good agreement is found. 



Abstract 



O ' The recent paper "There Are Infinitely Many Prime Twins" by R.F. Arenstorf [1] has 

2^ ■ raised a lot of excitement. The author claims to proved that: 

O; lim J2 log(p)log(p + 2) = C2 (1) 

p<N ,p,p+2 twins 

3 ! where the twin constant 

C2 = 2Y\(l- 7 — ^—^] = 1, 32032363169373914785562 .... (2) 

V; p>2^ (p-iyj 

H 

CZ \ Waiting for the formal approval of this result by mathematical community I have run 

the computer program to check the validity of (1). Even Hardy and Littlewood in their 

famous paper [2] have presented tables with numerical verification of their conjectures 

based on the existing that time data up to 9 000 000. Because the program I have 

written several years ago performs all operations on bits it was natural to store the data 

representing the actual arithmetical mean value of log(p) log(p + 2) at values of N forming 

the geometrical progression with the ratio 2, i.e. at A^ = 2^^, 2^^, . . . , 2^^, 2^°. The results 

are presented in the Table 1. 



TABLE I 



N 


l/iVEp<ivlog(p)log(p + 2) 


l/Nj:,<N^og{p)\og{p + 2)/C2 


222 ^ 4194304 


1.330875543 


1.00799191245 


2^^ = 8388608 


1.325154123 


1.00365856579 


224 ^ 16777216 


1.323501313 


1.00240674443 


225 = 33554432 


1.320938577 


1.00046575325 


226 = 67108864 


1.319310330 


0.99923253570 


227 = 134217728 


1.320265943 


0.99995630684 


228 = 268435456 


1.319095515 


0.99906983694 


229 = 536870912 


1.319679380 


0.99951205023 


230 ^ 1073741824 


1.320096901 


0.99982827612 


2^^ = 2147483648 


1.320047000 


0.99979048210 


232 ^ 4294967296 


1.320350510 


1.00002035747 


233 _ 8589934592 


1.320423613 


1.00007572498 


234 ^ 17179869184 


1.320490713 


1.00012654539 


235 _ 34359738368 


1.320290443 


0.99997486310 


236 _ 68719476736 


1.320309503 


0.99998929874 


237 ^ 137438953472 


1.320365187 


1.00003147324 


238 ^ 274877906944 


1.320351882 


1.00002139677 


239 ^ 549755813888 


1.320340769 


1.00001297956 


240 = 1099511627776 


1.320322532 


0.99999916675 



As it can be seen from above table there is no apparent dependence on A^. Indeed, 
trying to find heuristically the dependence on N we can argue that the probabihty to find 
twin pair around x is C2/log^(a;) and hence the mean expectation value of the product 
log(p) log{p + 2) for p and p + 2 on both sides of a: {p = x — l,p + 2 = x + l) does not 
depend on x and we have simply that 



J2 \ogip)\ogip + 2) = C2N 

p<N,p,p+2 twins 

It can be contrasted with the calculation of the Brun constant 

1 i\ n 1\ f\ 1 



(3) 



^^n3njn5+7Jnii+i^j+-<-- (^) 

The probability to find a pair of twins in the vicinity of x is 2C2/log {x\ so the 
expected value of the finite approximation to the Brun constant can be estimated as 
follows: 



B,{x) = B,(oo) 



1 ^ _ f"^ du 
-^Bo- 4co 



ptwinprimOxP "' ""' -^^ u\og\u) 



Bo 



4C2 

log(a;) ■ 



(5) 



It means that the plot of finite approximations B2{x) to the original Brun constant is a 
linear function of l/log(a;) [3] and from the partial sum B{x) calculated on the computer 



up to X the limiting value can be extrapolated by adding 4C2/log(a;): B = B{x) + 
4C2/log(x). To gain some idea what value of the limit can be the extrapolated from 
numbers in Table I the Figure 1 presents actual values of the mean value of log(p) log(p+2) 
plotted against 1/A^. Fitting the straight line to these points by least square method gives 
the intercept (what corresponds to A^ = oo) 1.3200385787619. In fact we see in Table I 
shortage of twins in the interval (2^^, 2^^) and in the next intervals some surplus of twins. 
Thus skipping the first 10 points and fitting straight line in the interval (2^^, 2^°) (in fact 
only two points are needed to determine straight line!) I got for the limiting value of the 
intercept 1.3203501777. 
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Figure 1: The plot of the actual mean values of log(p) log(p + 2) against 1/A^. Notice that 
on the X axis there is a logarithmic scale. 



